Textures of Spin-Orbit Coupled F = 2 Spinor Bose Einstein Condensates 
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We study the textures of F — 2 spinor Bose-Einstein condensates (BECs) with spin-orbit coupling 
(SOC) induced by a synthetic non-Abelian gauge field. On the basis of the analysis of the SOC 
energy and the numerical calculation of the Gross-Pitaevskii equation, we demonstrate that the 
textures originate from the helical modulation of the order parameter (OP) due to the SOC. In 
particular, the cyclic OP consists of two-dimensional lattice textures, such as the hexagonal lattice 
and the |-vortex lattice, commonly understandable as the two-dimensional network of the helical 
modulations. 
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Motivated by the recent successfulgeneration of gauge 
fields in neutral cold atoms 0, 0, much attention 
has beenpaid to Bose-Einstein condensates (BECs) in 
Abelian [ij and non-Abelian [H-Q gauge fields. The syn- 
thetic field technique has allowed access to the quantum 
Hall regime [5| and has realized nontrivial textures in 
spinor BECs 

It is well known in fcrmion systems such as topological 
insulators and superconductors that spin-orbit coupling 
(SOC) plays a crucial role. This is particularly true for 
fermionic superfluid 3 He because SOC due originally to 
the 3 He nuclear dipole-dipole interaction is decisive in de- 
termining the order parameter (OP) textures, although 
its force is so tiny that it modifies the Coop er p air struc- 
ture itself, as first pointed out by Leggett [1(1. Since in 
bosonic superfluids, namely the spinor BEC [ill. [l2|, no 
corresponding "natural" SOC exists, one may implement 
it by using a "synthetic" gauge potential to better under- 
stand and control the OP space in spinor BEC systems. 
In fact, we see that by introducing SOC, they exhibit a 
plethora of types of the OP textures in a controlled way, 
such as fractional vortices, two-dimensional periodic tex- 
tures, and skyrmions. 

The pscudospin states are the eigenenergy states of 
the laser-atom interaction, sufficiently isolated from the 
other states. The adiabatic motion of these states pro- 
duces a vector potential in real space EH • Several 
schemes to generate gauge fields are suggested. One of 
them is the A scheme 0,3! which is useful to realize a 
symmetric spin-orbit coupling term, such as the Rashba- 
type SOC. The different scheme is implemented by the 
NIST group (H, where the SOC consists of the equal con- 
tribution of the Rashba and Dresselhaus types. 

The gauge transformation cannot remove the spatially 
uniform non-Abelian vector potential. Thus, the BECs 
have a spatially modulated OP and textures due to finite 
momentum [6H8| . Nontrivial textures due to the Rashba- 
type SOC in pseudospin F = 1/2 and F = 1 BECs have 
been studied recently Q- The two possible textures are 
the one-dimensional plane and standing waves; in the 
former, the phase of the OP varies and in the latter, 



the amplitude oscillates along a favorable direction. In 
F=l BECs, the stable region of the plane and standing 
waves is equivalent to that of the ferromagnetic and polar 
phases, respectively [ri| . However, we emphasize here 
that an open question remains of how textural structures 
can emerge in F=2 spinor BECs, where the cyclic phase 
distinct from the manifold of the ferromagnetic and polar 
phase can be the magnetic ground state. 

The aim in this Rapid Communication is to clarify the 
role of the SOC in bosonic superfluids. We focus on F = 2 
spinor BECs, since they have magnetic ground states 
that include the cyclic phase. On the basis of the analysis 
of the SOC energy, we reveal that the plane (standing) 
wave realized in ferromagnetic (polar) phase can be in- 
terpreted as the rotation of the OP in pseudospin space, 
which simultaneously propagates along one direction in 
real space. In addition, skyrmions of the uniaxial polar 
OP emerge in the small size of the system. This analy- 
sis enables us to understand the emergence of nontrivial 
textural structures in the cyclic phase. We demonstrate 
that the hexagonal lattice and |-vortex lattice states are 
energetically competitive. We also calculate the phase 
diagram where the spin-spin interaction favors the cyclic 
phase. 

We consider a zero temperature F = 2 spinor BEC 
with a SOC. The single-particle Hamiltonian including 
the Rashba-likc SOC term is described as 
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where $f(r) = [i/>2, tpi, ipo, V'-ii "0-2] T is the OP vector 
m an F = 2 BEC, /i = -V 2 /2m + Vp 0t (r) consists of the 
kinetic energy and the trap potential term, and e± = x ± 
iy. The planar hexapod setup introduced by Juzeliiinas 
et al. [H realizes Rashba-like SOC with ,5 = 1. However, 
in an F = 2 spinor BEC system, the experimental setup 
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for realizing the precise Rashba-typc SOC form with f} = 
s/6/2 has not been found thus far. 

To reveal the role of the SOC term in Eq. (JXJ) , we first 
simplify the OP as 



#(r) = R(k ■ r,n R )*i(r ), 



(2) 



where R(9, ft) denotes the rotation matrix with the an- 
gle 9 about h in the pseudospin space; fin. and k are 
the rotation axis and the modulation vector in the x-y 
plane, respectively. ^\ is an arbitrary OP vector in the 
pseudospin space at a certain point r$. By substitut- 
ing Eq. ((2]) in Eq. ([1]), one can obtain the SOC energy 
density: 
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where yl = (1 + 76/3), S= ^/S/2 (3 + ^6(3), C= (3-^6/3), 
u; = exp[ifc • r], ^j = _R(7r/3, ein)$i, and em || x + y + z. 

The diagonal elements of Afg Q in Eq. ([3]) are found to 
be energetically dominant, compared to the off diagonal 
elements. This is because the oscillation terms with w™ 
(ti^O) hardly contribute to energetics when they are in- 
tegrated over r. The dominant elements are proportional 
to k ■ n R . (i) The SOC favors the situation where the ro- 
tation axis fiR corresponds to the modulation vector k. 
We call this spatial modulation the helical modulation. 
In addition, we define the angular resolved SOC energy 
density hAn(k, Vl^i) = hso(k, fin = k, ^i), which denotes 
the SOC energy density for with a given modulation 
vector k. (ii) The helical modulation of the OP is inde- 
pendent of the details of SOC term, such as the precise 
Rashba-type SOC ((3 = y/6/2, that is, C = 0) and the 
Rashba-like SOC ((3 = 1) implemented by the planar 
hexapod setup 

Here, we assume that the atoms in pseudospin states 
interact through the most symmetric interaction: 



Hint = \ j d 2 r [c n 2 + Cl S ■ S + c 2 \A 00 \ 2 ] , 



(4) 



where n = #+\&, S = ^crfy, and A oa = (2^2^-2-2^1^-1 + 
V'o' ! /'o)/v / 5 are the particle density, the spin density, and 
the singlet pair amplitude, respectively. The purpose of 
using this assumption is to concentrate on the symme- 
try breaking due to the SOC term. In addition, it is 




FIG. 1: (Color online) OP profiles in pseudospin space and 
energetically favored modulation vectors of ground states: (a) 
FM, (b) BP, and (c) UP. The color on the surface corresponds 
to the phase of OP. The UP-skyrmion texture is shown in (d). 



known that the minimization of H mt leads to four mag- 
netic ground states: Ferromagnetic (FM), biaxial polar 
(BP), uniaxial polar (UP), and cyclic (CY) phase 12 1. 

We numerically minimize the full Gross-Pitaevskii en- 
ergy functional Hq+H^ without any restriction by using 
the imaginary time evolution scheme in the presence of 
the cylindrical symmetric trap. The numerical results 
presented in this Rapid Communication are obtained by 
setting the parameter (3 =1 and all of the results are not 
changed qualitatively by the parameter (3. The numeri- 
cal calculation reveals that for the FM parameter region, 
as shown in Fig.[TJa), the resulting OP can be simplified 
to the one-dimensional helical modulation of the uniform 
FM OP along the y axis, that is, = R(ky, y)^FM 
with vj/ FM = i?- 1 ( 7 r/3,eiii)[l,0,0,0,0] T . The wave num- 
ber k depends on k. In the same way, the favorable 
OP for the BP is given by replacing $fm m ^bp = 
R- 1 (tt/3, e m ) [1, 0, 0, 0, 1] T , which is shown in Fig. [T|b). 
Using these simplified OPs, the SOC energy density re- 
duces to \h AR (ky, *fm(bp))I = h 2 nk(V6(3 + l)/4-\/2, in- 
dependent of the FM and BP. These one dimensional he- 
lical modulations of the FM and BP are consistent with 
those realized in F = l/2 and F = 1 systems Q ■ 

In contrast, the favorable OP in the UP state reduces 
to $(r) = R(k ■ r,fe)* UP with vj/ UP = [0, 0, 1, 0, 0] T , 
where the modulation vector k can be arbitrary in the 
x-y plane as shown in Fig. QJc). Because k is not unique, 
the most favorable modulation turns out to be k \\ r, 
resulting in the skyrmion texture in Fig. [Tf d) . The 
SOC energy density of the UP skyrmion is found to be 
\h AR (kr,^ijp)\ = 3h 2 Kk(V6(3 + l)/16y/2, which is less 
than that of the FM and BP. 

We note that the UP and BP are degenerate in the 
absence of SOC K = 0. In the presence of the Rashba-like 
SOC, the energetics of the UP skyrmion texture and BP 
standing wave is determined by the ratio of two length 
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scales 1/k and the Thomas- Fermi radius -Rtf, which 
characterize the modulation in real space and the size 
of the condensate, respectively. The stable region of the 
UP skyrmion is 1/k > Rtf > while the BP standing wave 
becomes stable for 1/k<Rtf- 

Next, we move to the textures in the CY parameter 
region with d > and c 2 > 0. The OP of the CY is 
characterized by a nodal structure as shown in inner bot- 
tom of Fig. which has eight point nodes atzbr±y±£, 
±x^y±z,zbx±y^fz, and ±x^ip^z and six antinodes. Their 
point nodes are connected with saddle lines. The inher- 
ent difference from those shown in Fig. [1] arises from the 
three-dimensional form of the nodal points. This three 
dimensionality gives rise to the textures with the two- 
dimensional modulation vectors as mentioned below. 

Before going to the numerical results, we discuss the 
role of the SOC term in Eq. ((I) for the cyclic phase. 
It is convenient to introduce two OP vectors = 
R{axccos[l/\/2,], x + y)^ CY = *p and *i = R(n/2,x + 
y)# C Y = *s- Here, * C y = [«A °> V \/2, 0, i/2] denotes 
the simple form of the cyclic OP shown in inner bottom 
in Fig. [2] and (v&s) describes the cyclic OP where one 
of point nodes (saddle points) points to the z axis. The 
two upper panels in Fig.[2]show the angular resolved SOC 
energy |/i A R(fc,*s)| and \h AR (k,$ P )\ for a given |fc|. In 
the case of $s ; the OP has four point nodes and two 
antinodes in the x-y plane. It is seen in the upper panels 
of Fig. [2] that for the case of |^AR,(fc, ^s)| becomes 
maximum (minimum) when k points in the nodal (antin- 
odal) direction. In contrast, the angularly resolved SOC 
energy density for h\R.(k, ^p) becomes cylindrically sym- 
metric in the x-y plane. The main panel of Fig. [5] shows 
£ar(#,¥>) = / d(j) k h AR (k^i) with = i?(i9,e v )* C Y 
and e v = — sin ipx + cos ipy. Here, it is found that the 
most stable situation in the cyclic region is = \I/p, 
where one of the point nodes points to the z direction. 
However, note that the SOC in Eq. (TTJ) also gives rise 
to the helical modulation of the OP, which rotates the 
direction of the point node from the z axis. 

Let us argue the stable textures in real space for 
ci/co>0 and C2/cq>0. Based on the numerical calcula- 
tion of the full CP equation including the SOC term, we 
find three kinds of textures: The uniform cyclic texture, 
the CY-UP hexagonal lattice texture shown in Fig. [31 
and the ^-vortex lattice texture shown in Fig. [4] 

The emergence of the CY-UP lattice texture is 
straightforwardly comprehensible, on the basis of the ar- 
gument of the SOC energy in the cyclic OP. As shown in 
Figs.EJa) and[3]Jb), this texture consists of the hexagonal 
unit cell in which the polar state is localized. Here, we 
look carefully into the texture inside the hexagonal unit 
cell. As seen in Fig.[3tc), the unit cell consists of the uni- 
axial polar OP at Uq enclosed by cyclic OPs at Pi to P 6 , 
which arc obtained from the numerical calculation of the 
GP equation. At all the points P\-P&, one of the point 




FIG. 2: (Color online) The main panel describes the SOC 
energy EAR(-d,ip) = J d(j>khAR{k, >5i). The bottom inner fig- 
ures show the cyclic OP V&cy and the angles § and tp, which 
denote the z axis of vPi. The two upper panels are the angu- 
lar resolved SOC energy \hAn(k, %)| for the cases of \&i = $s 
and tyi = typ. In these panels, the radii of the outer and inner 
lines indicate |/iAa(fe, *Pi)| and the deviation of \hAB.(k, 
from the value at the antinode, in the k x -k y plane. 




FIG. 3: (Color online) Profiles of the CY-UP lattice texture: 
(a) pseudospin density S(r), (b) singlet amplitude Aoo(r), 
and (c) the order parameter. The parameters are set to be 
k7?tf = 10.1 and ci/co = 02/00 —0.2. 



nodes always points to the z axis, which makes the SOC 
energy lower, as discussed in Fig.[5J In addition, at Pi the 
other three point nodes are located at (6, <j>) = (9 , — tt/3), 
(6*o,7t/3), and (0o,tt) with 60 = arccos(2v / 2/3). By prop- 
agating the helical modulation around one of the saddle 
points, the OP $p at Pi can be continuously transformed 
to that at P2, which orients the point node (9q, — 7r/3) at 
Pi to the z direction at P2. The point node (0q,tt/3) 
and (60 , 7r) at Pi is oriented to the z direction at Pq 
and P4. In the same way, the OP at Pi can propa- 
gate with the helical modulation along the circumference 
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FIG. 4: (Color online) Profiles of the |-vortex lattice texture: 
(a) pseudospin density S(r), (b) singlet amplitude Aoo(r), 
and (c) the order parameter. The parameters are set to be 
Ki?TF = 10.1, ci/co=0.2, and C2/co = 20. 




FIG. 5: (Color online) Phase diagram of textures spanned by 
coupling constant of SOC k and that of the spin-spin interac- 
tion C2/C0. We fix the parameter ci/co = 0.2. The stable phase 
for kRtf < 1 continuously changes from the UP skyrmion to 
uniform CY texture. In the kRtf > 1 region, three textures 
appear: The BP standing wave for the small C2/C0, the CY- 
UP lattice (hexagon), and i-vortex lattice (square). 



P2 —t P3 —> Pi . Note that the cyclic phase at Pi can 

have six modulation vectors which have isotropic energy 
gain in the SOC term, as shown in Fig. [5J The central 
region of the unit cell in Fig. [3Jc) is occupied by the UP 
OP to avoid frustration of the cyclic OP at Uq. Hence, 
in this sense this network of the helical modulation can 
be regarded as the close-packed hexagonal lattice of the 
cyclic OP whose point node points to the z axis. 

The other possible texture in the cyclic region is dis- 
played in Fig. IH which can be the ferromagnetic core at 
Fq enclosed by cyclic OPs. At A\ in Fig. EJc), the antin- 
ode points to the z axis, which cannot be the ground state 
of the SOC energy in Eq. (P). The OP at A x has fourfold 



symmetry in the x-y plane. The OP at A\ can be contin- 
uously transformed to that at A2 by the helical modula- 
tion along ~x — y with the rotation angle ir/2. After the 
helical modulation along the path A-y —tA^ — >■ ■ -A\ in the 
same way, the OP results in the shape where the antinode 
points to the z axis. However, the resulting OP differs 
from the original one and has the phase shift +2tt/3. In 
order to recover the single- value nature at A\ , the phase 
shift — 27r/3 is compensated by the i7(l) phase of the 
OP. Hence, this texture with fourfold symmetry can be 
regarded as the square lattice of the | vortices, in the 
sense that the U(l) phase of the cyclic OP continuously 
changes by 2ir/3 along the path A± — »• A2 — >• ■ ■ • — >■ A\. 
Here, the central region at Fq is occupied by the ferro- 
magnetic core because of frustration of the cyclic OP. It 
is expected that the ^ vortex also appears in the cyclic 
phase of F = 2 spinor BEC under rotation [l4j|, which 
behaves as the non-Abelian vortices. 

In Fig. [5j we summarize the phase diagram spanned 
by kRtv and C2/C0. For kRtf ;$ 1, the stable texture 
continuously changes from the UP skyrmion to the uni- 
form CY. For kRtf ^ 1, the standing wave of the BP 
appears in the small C2/C0 region. As C2/C0 increases, 
the hexagonal lattice of UP core and the ^-vortex lattice 
become energetically competitive. The former texture is 
the close-packed lattice of the CY OP with the point node 
pointing to the z axis. On the other hand, the FM core 
in the ^-vortex lattice relatively favors the Hi nt , which 
becomes stable in the large c-i / cq region. 

In summary, we have studied the stable textures of 
F = 2 spinor BEC with a spin-orbit interaction. Based 
on the detailed analysis of the spin-orbit interaction, we 
find that the ferromagnetic and biaxial polar OPs have a 
unique direction for the helical modulation, while that of 
the uniaxial polar OP is not unique. The nonuniqueness 
leads to the emergence of the skyrmion textures, which 
become stable in the parameter region where the polar 
phase is stable. Moreover, we have computed the phase 
diagram where the cyclic phase is the magnetic ground 
state. The phase diagram is covered by two-dimensional 
lattices, such as the CY-UP lattice and i-vortex lat- 
tice. These lattice structures are understandable as the 
two-dimensional network of the helical modulation of the 
cyclic order parameters. We emphasize that those novel 
textures arc not seen in rotating F = 2 spinor BEC [l5[ 
and appear inherently due to SOC. 

We should mention that the Hamiltonian for spinor 
BECs under synthetic gauge fields does not necessarily 
have the pseudospin rotation symmetry Q, which can 
be altered by the laser configuration and so on. Hence, 
the general form of the Hi nt is open for future study. 
Moreover, the experimental observation of the textures 
remains a challenging problem different from those with 
hyperfine spins [161 ] 

Note added: Recently, we became aware of a preprint 
by Xu et al. [l7| . which has some overlap with our results. 
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